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ABSTRACT 

We use stellar masses, surface photometry, strong lensing masses, and stellar velocity dispersions 
(cc/2) t° investigate empirical correlations for the definitive sample of 73 early-type galaxies (ETGs) 
that are strong gravitational lenses from the SLACS survey. The traditional correlations (Fundamental 
Plane [FP] and its projections) are consistent with those found for non-lens galaxies, supporting 
the thesis that SLACS lens galaxies are representative of massive ETGs (dimensional mass M<ji m = 
10 11 — 10 12 Mq). The addition of high-precision strong lensing estimates of the total mass allows us 
to gain further insights into their internal structure: i) the average slope of the total mass density 

profile (ptot oc r~ 7 ) is (7') = 2.078 ±0.027 with an intrinsic scatter of 0.16 ±0.02; ii) 7' correlates with 
effective radius (r c ) and central mass density, in the sense that denser galaxies have steeper profiles; 
iii) the dark matter fraction within r c /2 is a monotonically increasing function of galaxy mass and size 
(due to a mass-dependent central cold dark matter distribution or to baryonic dark matter - stellar 
remnants or low mass stars - if the IMF is non- universal and its normalization increases with mass); 
iv) the dimensional mass Maim = 5r c er 2 , 2 /G is proportional to the total (lensing) mass M ro /2, and 

both increase more rapidly than stellar mass M* (M* ex M°%); v) the Mass Plane (MP), obtained by 
replacing surface brightness with surface mass density in the FP, is found to be tighter and closer to 
the virial relation than the FP and the M*P, indicating that the scatter of those relations is dominated 
by stellar populations effects; vi) we construct the Fundamental Hyper-Plane by adding stellar masses 
to the MP and find the M* coefficient to be consistent with zero and no residual intrinsic scatter. 
Our results demonstrate that the dynamical structure of ETGs is not scale invariant and that it is 
fully specified by M r / 2 , r c , and <t c /2- Although the basic trends can be explained qualitatively in 
terms of varying star formation efficiency as a function of halo mass and as the result of dry and wet 
mergers, reproducing quantitatively the observed correlations and their tightness may be a significant 
challenge for galaxy formation models. 

Subject headings: galaxies: elliptical and lenticular, cD - galaxies: fundamental parameters - galaxies: 
structure - dark matter - gravitational lensing 



1. INTRODUCTION 

The hierarchical model for structure formation in the 
context of a cold dark matter cosmology (ACDM) has 
been tremendously successful at describing the larg e 
scale features of the Universe (e.g.. lKomatsu et aI1l2009l) . 
However, there are many important properties of the 
Universe at galactic and sub-galactic scales that have 
evaded a detailed understanding. For example, the sev- 
eral bi-modal classes of ga laxies (e.g., red/bl ue color, 
early/late- type morphology: iBalogh et al.ll2004|) . the so- 
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called downs izing of star forma tion (iCowie et al.l Il996t 
iCooper et aHl2006t iBundv et al.l l2006[). the tight empiri- 
cal correlations between the observed properties of early- 
type galaxies (Fab er fc Jackson! [1976 : Kormendy 1977; 
iDressler et~aT1 11987: Diorgovski & Davislll987t) . correla- 
tions (or the l ack of correlations) between these an d the 
local density (|Dresslerl U980t ICooper et aTll2006D . and 
the absence of lo cal analogs to extremely compact high- 
redshift galaxies (IDaddi et alj |2005t iTruiillo etalll2006t 
Ivan Dokkum et al.ll2008ft . 

One important step in understanding these phenom- 
ena is to explore the relationship between the baryonic 
matter that dominates astrophysical observables and the 
dark matter that is postulated in the ACDM model. High 
precision measurements at galaxy scales are essential to 
test whether apparent inconsistencies between the ob- 
served universe and dark matter only cosmological simu- 
lations can be reconciled by an improved understanding 
of the detailed physical mechanisms governing baryons 
and their interaction with DM, or whether a rethinking 
of the CDM paradigm might be necessary. Within this 
context, the origin of early- type galaxies (ETGs) is cur- 
rently a point of discord between observation and theory, 
and therefore its investigation carries enormous potential 
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for discovery. Although their formation via merging of 
spiral disks is one of the assumptions of the standard 
paradigm, it remains to be seen whether this can work 
in detail. 

The tight empirical correlations between the observed 
properties of ETGs are a powerful phenomenological 
tool to relate baryonic and dark matter. Among these, 
the correlation between size, surface brightness and 
stellar velocity dispersion known as the Fundam ental 
Plane ([Faber et alJll987tlDiorgovski fc Davislll987L here- 
after FP), and i ts mass counterpart the Stellar Mass 
Plane (M*P; e.g.. lHvde fe Bernar di 2009b], where surface 
brightness is replaced by stellar mass) have provided two 
key insights. First, the correlations are "tilted" in the 
sense that they cannot be explained by assuming that 
ETGs are a self-similar family, they obey th e virial theo- 
rem, and have a constant mass-to- light ratio (jFaber et al.1 
Il987t iCiotti et al.l Fl996): mounting evidence, including 
the alignment of the Mass Plane (MP; IBolton et al.ll2007ft 
with the virial plane, suggests that the dominant cause of 
the ti lt is a mass-dependent central dark matt er fraction 
(e.g., iBolton et al1l2008bt iTortora et aT1l2009h. Second 
the correlations are remarkably tight (LJo rgcnse net all 
fl99fil:IHvde fc BernardillMMlGraves fc Faberll2010h im- 
plying that for a given size and velocity dispersion there 
is very little scatter in the star formation histories of 
ETGs. Previous studies have f ound a slight mis alignment 
between the FP and M*P (jHvde fc Bernardi! 12009b') 
which likely results from constructing the FP with galax- 
ies that span a variety of ages, since age correlates with 
both stellar mass (more massive galaxies form earlier) 
and the stellar mass-to-light ratio (older stellar popula- 
tions have larger stellar M/L for a given stellar mass). 

Additional progress in understanding ETGs has been 
made by attempting to infer the separate luminous 
and dark componen ts in the central regions of galaxies. 
ITortora et al.l ([2009ft use dynamical masses determined 
from central velocity dispersion measurements with stel- 
lar masses inferred from stellar populations synthesis 
(SPS) models to obtain the central dark matter frac- 
tion; they find a clear dependence on total mass in the 
sense that more massive galaxies have lar ger dark matter 
fractions (also see lNapolitano et"aT1l2010ft . although they 
assume a mass-traces light or singular isothermal sphere 
mass distribution. The assumption of an isothermal cen- 
tral mass profile seems to be robust, as illustrate d by 
modeling of X-ray data ([Humphrey fc Buotel 12009ft and 
stron g lensing and stellar kine matics (iKoopmans et al.l 
[200l 12001 iBarnabe et aT][2009HBarnabe et al.ll2010ft . al- 
though the origin of this isothermality remains unclear. 

We use a sample of early- type gravitational lenses from 
the Sloan Lens ACS Surv ey f SLACS: iBoTton et aLlEJOOl 
l2008at lAuger et al.ll2009ft to investigate the scaling rela- 
tions of ETGs, and in particular we look at the relation- 
ships between total and stellar mass with respect to the 
other structural parameters of lenses. We exploit multi- 
ple lines of observation - including strong lensing, stel- 
lar dynamics, and SPS models - to distinguish between 
the luminous and dark mass, and we construct the FP, 
M*P, and MP for the SLACS ETGs. The paper is or- 
ganized as follows. In Section 2 we present the SLACS 
dataset, summarize observables listed in previous papers 
and used for this analysis, and present new quantities, 
such as improved estimates for the total mass within half 



of the effective radius and the central slope of the total 
mass density profile. Section 3 describes bivariate empir- 
ical correlations derived for the SLACS sample, starting 
from traditional ones including non-lensing observables 
and concluding with those including total mass as derived 
from strong lensing. Section 4 describes correlations in 
higher dimensions, including the FP, M»P, MP and the 
newly derived Fundamental Hyper-Plane. Section 5 dis- 
cusses our results and Section 6 gives a brief summary. A 
standard cosmological model with Q m = 0.3, il\ = 0.7, 
and h = 0.7 is assumed throughout. 

2. THE SLACS SURVEY DATA: OBSERVABLES AND 
DERIVED QUANTITIES 

There are 85 confirmed (grade 'A') strong gravitational 
lenses that have been disc overed by SLACS ([Bolton et al.l 
l2008aH Auger et alJl2009l ). These include 73 galaxies with 
E or SO morphologies which we take to be the ETGs from 
SLACS and which we focus on in this paper. A wealth 
of data exists for each of these galaxies, including high- 
resolution multi-band optical and near-infrared Hubble 
Space Telescope (HST) imaging and fiber-based optical 
spectroscopy from the Sloan Digital Sky Survey (SDSS; 
lYork et all 120001) . These data are used to infer several 
fundamental properties of the lenses, which we briefly 
detail below and list in Table [1] Some of the observables 
and derived quantities have been given in previous papers 
of the SLACS series, and will not be repeated in Table [TJ 
for conciseness. The present compilation is based on the 
most up-to-date data and calibrations and supercedes the 
information presented in previous papers. 

2.1. Stellar Masses, Luminosities, and Effective Radii 

The high-resolution multi-band HST imag ing is used 
to in fer stellar masses M* for each system ([Auger et al.l 
2009:) using SPS models and assuming either a Chabrier 
or Salpeter IMF. Four systems only have one band of 
HST imaging and we exclude these from our analysis 
when stellar masses are required. The SPS models pro- 
vide robust synthetic photometry (sometimes referred to 
as fc-corrections, or k-color corrections) that is approxi- 
mately insensitive to the assumed IMF and these models 
therefore yield accurate estimates for the B- and F-band 
rest-frame luminosity. The same models also allow us to 
compute rest frame luminosities of each galaxy passively 
evolved to z = in a self consistent manner; B- and 
V- band luminosities at the redshift of the l ens and cor- 
recte d to z = can be found in the paper by ([Auger et al.l 
[2009ft . 

Furthermore, we use the HST imaging to determine 
the effective radius in each band, and we employ a lin- 
ear model of effective radius as a function of wavelength 
to infer the rest-frame F-band effective radius. We as- 
sume r e ,x = a* X + b where a and b are determined from 
a fit to the observations of r c in each filter with rest- 
frame wavelength given by A c /(1 + z) where A c is the 
filter central wavelength (e.g., iTreu et aLl 12001ft . Then 
r e,5500 is the effe ctive radius us e d thro ughout this paper. 
As discussed bv IBolton et al.l ([2008bf ). our assumption 
of Ide Vaucouleursl (| 1948ft profiles is valid in the lumi- 
nosity range covered by the SLA CS sample. The sys- 
tematic trends in Sersic index n (Sersic 1968) for such 
high-luminosity galaxies are dominated by the intrinsic 
scatter in the correlation. As expected, by fitting the 
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SLACS lenses with Sersic models, we find no correla- 
tion between n and any of the global galactic quantities. 
Likewise the adoption of Sersic profiles does not change 
any of the trends presented here. Therefore, we limit our 
analysis to the simpler and better constrained de Vau- 
couleurs models. This does not affect our interpretation 
of the tilt of the FP and other key structural parameters, 
as the effect of varying n is negligible in the probed mass 
range (jNipoti et al.ll2008h . 

2.2. Stellar Velocity Dispersions and Dynamical Masses 

The SDSS spectroscopy provides an estimate of the 
luminosity- weighted stellar velocity dispersion within the 
3"-diameter aperture of the SDSS fiber s which we refer 
to as ordss- We use the prescription of Uorgensen et all 
(f 995) to infer the velocity dispersion within half of the 
effective radius, cr e /2, for our analysis of the scaling rela- 
tions and parameter planes but use the aperture velocity 
dispersion in our analysis of the central mass profile (Sec- 
tion mni). 

We combine effective radii and stellar velocity disper- 
sions to construct a dimensional mass, defined as 



5tv^/2 
G : 



(1) 



where the number 5 is the canonical choice for the virial 
coefficient of dynam i cal masses of massive ETGs (e.g., 
iBernardi et"alll2003t ICappellari et "all 12009) . We note 
that the di mensional mass is n ot actually the dynamical 
mass (e.g-. IBolton et al]|2008b|) . but we choose this form 
to simplify comparison with dynamical masses. 

2.3. Einstein Radii, Mass Density Profile Slopes, and 
Strong Lensing Masses 

Singular isothermal ellipsoid (SIE) lens model fits to 
the HST data have b een used to derive Einstein radii 
for e ach lens systems (jBolton et al.l l2008at lAuger et al.l 
2009). The lens and source redshifts are known from 
the SDSS spectroscopy, and we can therefore infer the 
SIE velocity dispersion, which we refer to as osie (e.g., 
iBolton eta l. 2008a). Our SIE mass models also robustly 
constrain the total projected mass within the Einstei n 
radii to a precision of a few percent (jBolton et al.| [2008a). 

We use the available information to constrain power- 
law total mass distribution models for each lensing 

galaxy. The ma s s d istributi ons are defined as in. 

iTreu fc Koopmansl 1)20041 ) and IKoopmans et all (|2006l 
2009), with p oc r~ 7 , and these models are constrained 
using the mass within the Einstein radius, the SDSS 
stellar velocity dispersion, and de Vaucouleurs fits to 
the stellar light distribution. Deta ils of how t he fits 
are impl emented can be found in [Suvu et al.l (2009) 
(also see IKoopmans et aTll2006ll2009( ). although for this 
analysis w e only use our baseline models, character- 
ized by a iHernquisH ([1990) model for the stellar dis- 
tribution and no anisotropy of the stellar orbits. The 
isotro py assumption is consist e nt with c omplementary 
(e.g., ITreu fc Koopmansl |2004t IKoopmans et al.l t "2009) 
and more detailed (e.g.. iBarnabe et al.ll2009l ) investiga- 
tions of the SLACS lenses. Note, however, that while 
mild radial anisotropy would lead to a slightly sha l- 
lower inferred mass slope (e.g., IKoopmans et all 120091) . 
it would cause the inference on the mass within a central 



aperture to be over- or under-estimated, depending on 
whether the aperture is smaller or larger than the Ein- 
stein radius. We list the updated mass slopes 7' for all 
lenses with robust kinematic and lensing data in Table 
1; this extends and completes our analysis of power- law 
mass models for t he SLACS lenses based upon the SDSS 
spectroscopy (e.g.. IKoopmans et al.ll2006ll2009l ). We use 
these power-law models to infer the total projected mass 
within half of the effective radius (denoted M re / 2 ; this 
radius is chosen because it is well-matched to the typical 
Einstein radius and therefore leads to the smallest errors 
from extrapolating the power-law mass model) which is 
used in our analysis of the scaling relations and parame- 
ter planes of the lens galaxies. 

2.4. Mass to Light Ratios and Dark Matter Fractions 

We construct three different estimators of the central 
mass to light ratio, by computing the ratio between total 
(lensing) mass, dimensional and stellar mass and lumi- 
nosity within r c /2. The three ratios are referred to as 
the total mass-to-light ratio (or M r /2/L), the dimen- 
sional mass-to-light ratio (Mdim/L) and the stellar mass- 
to- light ratio (M*/L). When needed, we use the symbol 
M/L to refer to the three mass-to-light ratios collectively. 

The mass-to-light ratio is relevant for understanding 
observations but simulations are more readily understood 
in the context of dark (or stellar) mass frac tion. We use 
the st ellar masses derived from SPS models (|Auger et al.l 
l2009f ) in conjunction with the total mass within half of 
the effective radius determined from lensing and dynam- 
ics to infer the projected (i.e. within a cylinder) dark 
matter fraction, /dm = 1 — M*/M ro / 2 where M* has been 
scaled to the stellar mass within r c /2 (this is 32% of the 
total stellar mass for a de Vaucouleurs distribution) and 
any gas is effectively treated as dark matter. We note 
that the projected DM fraction is always larger than the 
(three-dimensional) DM fraction within a sphere of equal 
radius because of the contribution of the outer parts of 
the halo to the projected quantity. We adopt in this 
paper the projected (two-dimensional) dark matter frac- 
tion because it is the most robustly determined quantity 
and the one closest to the observables. However, the 
interested reader can derive the three-dimensional DM 
fraction from the projected fraction using the effective 
radii and mass density profile slopes given in Table [T] 

3. BIVARIATE CORRELATIONS 

Bivariate empirical c orrelations between pro perties of 
ET Gs, e.g., the L - a (|Faber fc Jackson! 1 19761) and L - 
r e (|Kormendvl 119771 ) relations, are extremely useful for 
a number of practical applications, including modeling 
of complex data and simulations of mock catalogs. Al- 
though ETGs are known to be at least a two parame- 
ter family, the bivariate correlations encode critical in- 
formation about the distribution of ETGs in higher- 
dimensional parameter spaces and therefore can be used 
to provide additional tests of theoretical models. For ex- 
ample, although dry mergers generally move ETGs inside 
the Fundamental Plane, they also tend to move galaxies 
away from its two dimensional pr ojections ([Nipoti et al.l 
12003 IBovlan-Kolchin et all 12009) . 

We fit several bivariate empirical correlations to the 
SLACS data described in Section 2. We have determined 
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TABLE 1 

Mass and Structural Parameters for SLACS Early- type Lenses 





r c,V 

(kpc) 


°"e/2 

(km s~ x ) 


log[ M, 


:,/M ] 


log[ M re/2 /M ] 


/dm 




Lens Name 


Chab 


Salp 


Chab 


Salp 


V 


SDSSJ0029-0055f 


8.36 


231 ± 18 


11.33 ±0.13 


11.58 ±0.13 


11.13 ±0.03 


0.47 ±0.15 


0.04 ±0.28 


2.38 ±0.23 


SDSSJ0037-0942 


7.44 


282 ± 10 


11.48 ±0.06 


11.73 ±0.06 


11.36 ±0.02 


0.57 ± 0.07 


0.25 ±0.11 


2.14 ±0.07 


SDSSJ0044+0113 


6.12 


267 ± 13 


11.23 ±0.09 


11.47 ±0.09 


11.13 ±0.07 


0.59 ±0.10 


0.29 ±0.18 


2.31 ±0.24 


SDSSJ0216-0813 


13.19 


334 ± 23 


11.79 ±0.07 


12.03 ±0.07 


11.76 ±0.02 


0.65 ± 0.06 


0.40 ±0.11 


2.09 ±0.20 


SDSSJ0252+0039 


5.68 


170 ± 12 


11.21 ±0.13 


11.46 ±0.13 


10.97 ±0.03 


0.42 ±0.18 


-0.03 ± 0.32 


1.57 ±0.12 


SDSSJ0330-0020 


6.23 


220 ± 21 


11.35 ±0.09 


11.58 ±0.09 


11.14 ±0.05 


0.46 ±0.12 


0.09 ±0.21 


1.91 ±0.18 


SDSSJ0728+3835 


5.86 


219 ±11 


11.44 ±0.12 


11.69 ±0.12 


11.12 ±0.02 


0.31 ±0.19 


-0.23 ±0.37 


1.86 ±0.10 


SDSSJ0737+3216f 


14.10 


338 ± 16 


11.72 ±0.07 


11.96 ±0.07 


11.52 ±0.03 


0.49 ± 0.09 


0.10 ±0.16 


2.68 ±0.12 


SDSSJ0819+4534f 


7.63 


227 ± 15 


11.15 ±0.08 


11.40 ±0.08 


11.16 ±0.04 


0.68 ± 0.07 


0.43 ±0.12 


2.16 ±0.25 


SDSSJ0822+2652 


7.64 


263 ± 15 


11.43 ±0.13 


11.69 ±0.13 


11.32 ±0.02 


0.56 ±0.13 


0.21 ±0.23 


2.12 ±0.14 


SDSSJ0912+0029 


11.69 


322 ± 12 


11.71 ±0.07 


11.96 ±0.07 


11.71 ±0.02 


0.67 ± 0.06 


0.43 ± 0.09 


1.98 ±0.09 


SDSSJ0935-0003f 


20.09 


391 ± 35 


11.72 ±0.07 


11.96 ±0.07 


11.81 ±0.14 


0.73 ± 0.09 


0.53 ±0.15 


2.44 ± 0.36 


SDSSJ0936+0913 


7.00 


246 ± 11 


11.43 ±0.12 


11.68 ±0.12 


11.18 ±0.02 


0.40 ±0.17 


-0.07 ±0.30 


2.24 ±0.12 


SDSSJ0946+1006 


9.08 


265 ± 21 


11.34 ±0.12 


11.59 ±0.12 


11.43 ±0.02 


0.73 ± 0.07 


0.51 ±0.14 


2.01 ±0.18 


SDSSJ0956+5100 


8.58 


338 ± 15 


11.56 ±0.09 


11.81 ±0.08 


11.52 ±0.02 


0.64 ± 0.07 


0.36 ±0.13 


2.30 ± 0.09 


SDSSJ0959+0410 


3.34 


203 ± 13 


10.91 ±0.07 


11.15 ±0.06 


10.76 ± 0.02 


0.54 ± 0.08 


0.20 ±0.13 


2.05 ±0.15 


SDSSJ0959+4416 


7.27 


248 ± 19 


11.47 ±0.12 


11.72 ±0.12 


11.23 ±0.02 


0.42 ±0.17 


-0.03 ± 0.30 


2.14 ±0.21 


SDSSJ1016+3859 


4.38 


254 ± 13 


11.23 ±0.12 


11.48 ±0.12 


11.04 ±0.02 


0.48 ±0.15 


0.08 ±0.26 


2.19 ±0.11 


SDSSJ1020+1122 


6.23 


290 ± 18 


11.54 ±0.12 


11.80 ±0.12 


11.34 ±0.03 


0.46 ±0.16 


0.04 ± 0.28 


2.08 ±0.12 


SDSSJ1023+4230 


5.97 


247 ± 15 


11.33 ±0.12 


11.57 ±0.12 


11.19 ±0.03 


0.55 ±0.13 


0.20 ±0.23 


2.01 ±0.11 


SDSSJ1029+0420 


3.02 


215 ± 9 


11.04 ±0.12 


11.29 ±0.11 


10.71 ±0.02 


0.28 ±0.19 


-0.28 ±0.34 


2.28 ±0.10 


SDSSJ1106+5228 


3.61 


266 ± 9 


11.13 ±0.06 


11.37 ±0.06 


10.91 ±0.02 


0.47 ±0.08 


0.07 ±0.13 


2.40 ±0.07 


SDSSJ1112+0826 


6.48 


328 ± 20 


11.48 ±0.09 


11.73 ±0.08 


11.43 ±0.03 


0.63 ± 0.08 


0.34 ±0.13 


2.21 ±0.10 


SDSSJ1134+6027 


5.26 


243 ± 11 


11.26 ±0.12 


11.51 ±0.12 


11.06 ±0.02 


0.48 ±0.14 


0.06 ±0.25 


2.20 ±0.11 


SDSSJ1142+1001 


6.99 


225 ± 22 


11.30 ±0.08 


11.55 ±0.08 


11.22 ±0.02 


0.60 ± 0.08 


0.31 ±0.13 


1.90 ±0.23 


SDSSJ1143-0144 


10.25 


263 ± 5 


11.36 ±0.09 


11.60 ±0.09 


11.50 ±0.03 


0.77 ±0.05 


0.59 ±0.09 


1.92 ±0.06 


SDSSJ1153+4612 


4.00 


233 ± 15 


11.08 ±0.13 


11.33 ±0.13 


10.90 ±0.03 


0.50 ±0.14 


0.10 ±0.26 


2.28 ±0.13 


SDSSJ1204+0358 


4.59 


274 ± 17 


11.20 ±0.07 


11.45 ±0.06 


11.09 ±0.03 


0.58 ±0.07 


0.26 ±0.12 


2.29 ±0.11 


SDSSJ1205+4910 


9.04 


282 ± 13 


11.48 ±0.06 


11.72 ±0.06 


11.42 ±0.02 


0.63 ± 0.06 


0.36 ±0.10 


2.16 ±0.12 


SDSSJ1213+6708f 


6.51 


292 ± 11 


11.24 ±0.10 


11.49 ±0.09 


11.17 ±0.02 


0.61 ± 0.08 


0.31 ±0.14 


2.49 ± 0.08 


SDSSJ1218+0830 


7.62 


218 ±10 


11.35 ±0.08 


11.59 ±0.08 


11.26 ±0.02 


0.60 ± 0.08 


0.30 ±0.13 


1.82 ±0.11 


SDSSJ1250+0523 


6.88 


256 ± 14 


11.53 ±0.07 


11.77 ±0.07 


11.20 ±0.02 


0.31 ±0.11 


-0.22 ±0.20 


2.30 ±0.12 


SDSSJ1306+0600 


6.12 


241 ± 17 


11.19 ±0.08 


11.43 ±0.08 


11.22 ±0.02 


0.70 ± 0.06 


0.47 ±0.10 


1.89 ±0.14 


SDSSJ1313+4615 


6.51 


266 ± 18 


11.33 ±0.09 


11.58 ±0.08 


11.27 ±0.02 


0.63 ± 0.08 


0.34 ±0.12 


2.06 ±0.14 


SDSSJ1318-0313 


14.05 


211 ±18 


11.43 ±0.09 


11.67 ±0.09 


11.60 ±0.02 


0.78 ± 0.04 


0.61 ±0.08 


1.64 ±0.15 


SDSSJ1330-0148 


1.39 


194 ± 9 


10.43 ± 0.06 


10.67 ±0.06 


10.31 ±0.03 


0.57 ±0.07 


0.26 ±0.12 


2.25 ±0.10 


SDSSJ1402+6321 


8.92 


268 ± 17 


11.55 ±0.07 


11.79 ±0.06 


11.45 ±0.02 


0.59 ±0.07 


0.28 ±0.11 


1.97 ±0.14 


SDSSJ1403+0006 


5.10 


218 ±17 


11.20 ±0.08 


11.44 ±0.08 


10.97 ±0.02 


0.45 ±0.11 


0.04 ±0.19 


2.14 ±0.23 


SDSSJ1416+5136 


5.92 


248 ± 25 


11.40 ±0.08 


11.64 ±0.08 


11.22 ±0.05 


0.50 ±0.11 


0.12 ±0.20 


1.90 ±0.16 


SDSSJ1420+6019 


2.56 


208 ± 4 


10.93 ±0.06 


11.17 ±0.06 


10.59 ±0.02 


0.30 ±0.10 


-0.23 ±0.18 


2.28 ±0.07 


SDSSJ1430+4105 


10.41 


325 ± 32 


11.68 ±0.12 


11.93 ±0.11 


11.64 ±0.02 


0.63 ±0.10 


0.35 ±0.18 


2.06 ±0.18 


SDSSJ1436-0000 


10.34 


226 ± 17 


11.45 ±0.08 


11.69 ±0.09 


11.39 ±0.02 


0.63 ± 0.07 


0.35 ±0.13 


1.88 ±0.19 


SDSSJ1443+0304 


2.64 


218 ±11 


10.87 ±0.06 


11.12 ±0.06 


10.67 ±0.02 


0.49 ± 0.08 


0.09 ±0.14 


2.31 ±0.12 


SDSSJ1451-0239 


5.83 


224 ± 14 


11.17 ±0.07 


11.39 ±0.06 


11.00 ±0.03 


0.52 ± 0.08 


0.19 ±0.13 


2.24 ±0.19 


SDSSJ1525+3327 


14.13 


265 ± 26 


11.78 ±0.09 


12.02 ±0.09 


11.72 ±0.02 


0.63 ± 0.08 


0.35 ±0.14 


1.77 ±0.20 


SDSSJ1531-0105 


7.54 


280 ± 12 


11.43 ±0.09 


11.68 ±0.09 


11.35 ±0.02 


0.61 ± 0.08 


0.30 ±0.14 


2.13 ±0.08 


SDSSJ1538+5817 


3.63 


194 ± 12 


11.03 ±0.08 


11.28 ±0.08 


10.80 ±0.03 


0.45 ±0.10 


0.02 ±0.18 


1.90 ±0.14 


SDSSJ1614+4522f 


8.18 


183 ± 13 


11.21 ±0.13 


11.47 ±0.12 


11.07 ±0.06 


0.54 ±0.16 


0.15 ±0.27 


2.00 ±0.29 


SDSSJ1621+3931 


8.85 


239 ± 20 


11.45 ±0.06 


11.70 ±0.07 


11.41 ±0.02 


0.64 ± 0.06 


0.36 ±0.10 


1.80 ±0.16 


SDSSJ1627-0053 


6.87 


295 ± 14 


11.45 ±0.09 


11.70 ±0.09 


11.30 ±0.02 


0.54 ±0.10 


0.19 ±0.17 


2.33 ±0.10 


SDSSJ1630+4520 


7.80 


281 ± 16 


11.61 ±0.07 


11.86 ±0.07 


11.43 ±0.03 


0.50 ±0.09 


0.13 ±0.15 


1.97 ±0.09 


SDSSJ1636+4707 


5.93 


237 ± 15 


11.38 ±0.08 


11.63 ±0.08 


11.13 ±0.02 


0.42 ±0.12 


-0.04 ±0.21 


2.09 ±0.14 


SDSSJ1644+2625 


4.49 


234 ± 12 


11.18 ±0.09 


11.43 ±0.08 


11.00 ±0.02 


0.50 ±0.10 


0.12 ±0.16 


2.09 ±0.10 


SDSSJ1719+2939 


4.32 


295 ± 15 


11.22 ±0.08 


11.46 ±0.08 


11.12 ±0.02 


0.59 ±0.08 


0.28 ±0.14 


2.36 ±0.09 


SDSSJ2238-0754 


5.78 


200 ± 11 


11.20 ±0.06 


11.45 ±0.06 


11.07 ±0.02 


0.57 ±0.07 


0.23 ±0.12 


1.79 ± 0.12 


SDSSJ2300+0022 


6.88 


284 ± 17 


11.40 ±0.07 


11.65 ±0.07 


11.36 ±0.02 


0.64 ± 0.06 


0.37 ±0.10 


2.06 ±0.13 


SDSSJ2303+1422 


9.30 


253 ± 16 


11.47 ±0.06 


11.71 ±0.06 


11.46 ±0.02 


0.67 ±0.05 


0.42 ± 0.09 


1.86 ±0.13 


SDSSJ2321-0939 


6.52 


246 ± 8 


11.35 ±0.08 


11.60 ±0.08 


11.20 ±0.02 


0.54 ± 0.09 


0.19 ±0.16 


2.02 ±0.08 


SDSSJ2341+0000 


10.59 


206 ± 13 


11.48 ±0.08 


11.73 ±0.08 


11.45 ±0.02 


0.65 ± 0.07 


0.38 ±0.12 


1.62 ±0.12 



Note. — Columns are: 1) Lens name; 2) effective radius, corrected to the rcst-framc V-band, in kpc; 3) central velocity dispersion within 
half of the effective radius; 4-5) total stellar mass assuming a Chabricr (col. 4) or Salpeter (col 5) IMF; 6) total mass within half of the 
effective radius, as determined from our power-law mass distribution models; 7-8) dark matter fraction within half of the effective radius for 
Chabricr and Salpeter IMFs; and 9) slope of the power-law mass distribution, p oc r~ ~ f . 
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linear fits to each correlation (these correlations are fre- 
quently between the logarithm of physical quantities and 
a linear fit therefore represents a power-law model) that 
account for the errors in both the dependent and inde- 
pendent variables as well as covariance between the mea- 
surement errors, and we also explicitly determine the in- 
trinsic scatter. We use a Pyt hon impleme ntation of the 
fitting technique proposed by [Kelly (2007), which uses a 
Bayesian framework to avoid biases introduced by inap- 
propriate choices for the prior distributions of the inde- 
pendent variables; we have found that this is particularly 
important when the errors on the independent variables 
are significantly larger than the errors on the dependent 
variables. 

We have restricted our analysis of these correla- 
tions to the early-type lenses, which we assume share 
similar formation and e volution histories. However, 
Uiang fc Kochanekl (|2007| ) suggest that the SLACS lenses 
are either not an homologous population or the stellar ve- 
locity dispersions ha ve significantly under-est imated sys- 
tematic errors (e.g.. iHvde k, Bernard i 2009a, who show 
that different velocity dispersion codes produce system- 
atically different results). We find that 6 galaxies in our 
sample are significant outliers of the hyper-plane rela- 
tion between size, velocity dispersion, stellar mass, and 
total mass given by Equation [3] (these systems are indi- 
cated in Table [IJ . We have investigated the nature of 
these outliers but find that they generally do not stand 
out from the other lensing galaxies; none of these dis- 
crepant objects have disky structure, nor do they have 
anomalous spectral features. Furthermore, we have de- 
termined the stellar velocity dispersions for each object 
using three independent codes but the codes do not find 
a significant difference between these six objects and the 
others. Nevertheless, we use an abundance of caution 
and exclude these six objects from the fits that include 
intrinsic scatter, as these objects tend to dominate those 
relations. 

We begin with "traditional" correlations between non- 
lensing observables in § 13.11 The main purpose of this 
section is to compare our inferred correlations with those 
inferred from samples of non-lens galaxies to test the hy- 
pothesis that SLACS lens galaxies are representative of 
the overall population of massive ETC Previous SLACS 
papers have investigated this issue based on a number 
of tests and have found no evidence for any difference 
between the SLA CS lenses and ET Gs with similar ve- 
locity dispersions dTreu et al.|[20 06: Bolt on et al.l l2008a; 
iTreuet aLlfeOOft I Auger et al.ll2009D . This study updates 
and extends some of those tests by considering the larger 
sample and including in the analysis relations based on 
stellar mass. 

We introduce lensing observables in § 13.21 to investi- 
gate correlations between stellar, total, and dimensional 
masses, as a means to constrain the virial coefficient, the 
initial mass function, and dark matter content of ETGs. 
In § 13.31 we study the distribution of slopes of the to- 
tal mass density profile 7', extending the analysis previ- 
ously published by lKoopmans et all (|2009f ). The goal of 
this analysis is twofold. From a galaxy formation point 
of view, the distribution of total mass density profiles 
constrains the relative distribution of baryons and dark 
matter and therefore constrains quantities such as the 
star formation efficiency. From the point of view of grav- 
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Fig. 1.— The a, 



/2 



M* relation for SLACS lenses. The solid black 



line is a linear fit to the relation (including scatter), the dotted 
black lines indicate the intrinsic scatter, the gray band indicates the 
quadrature-sum of the scatter and the u ncertainty on th e linear fit, 
the red dotted line is the linear fit from Hvdc & Bernardi (2009a), 
and the red dashed line is their quadratic fit. We find that the 
SLACS cr e /2-M, relation is consistent with the SDSS relations, 
although our formal fit is shallower due to our explicit treatment 
of intrinsic scatter. 

itational lensing studies, the distribution of 7' is an es- 
sential piece of information for inferences regarding, e.g., 
cosmological parameter s from grayitatiqnal time d elays 
and lens statistics (e.g.. lSuvu et al. 2009; Dobkc & King 
[2TM[Qguri|[2007l) . 

Finally, in § 13.41 we examine variations in central mass- 
to- light ratio and dark matter fraction with galaxy global 
properties again as a means to investigate the initial mass 
function and dark matter content of ETGs. 

3.1. Traditional (non-lensing) correlations 

The cr e / 2 -M* relation is shown in Figure [1] for a 
Chabrier IMF. The relation is slightly shallower than for 
SDSS galaxies in general (a = 0.18 ±0.03, 6 = 2.34 ±0.01 
with a e / 2 in units of km s _1 and M* in units of 1O 11 M0), 
largely due to explicitly including intri nsic scatter in the 
relation (Table [5]). Nevertheless, the IHvde fc Bernardi 
(2009a) relations are acceptable fits to the SLACS data; 
this is another indication that SLACS lenses constitute 
a velocity dispersion-selected subsample of the g eneral 
population of massive ETGs ([Bolton et a l. 2008a). 

We examine the r e -M* relation in Figure [2] We find 
that the SLACS lenses have a somewhat steeper relation 
(a = 0.81 ± 0.05, b = 0.53 ± 0.02 for a Chabrier IMF with 
r c , measured in kpc and M* in units of 10 11 Mqj see 
Table Et than non-lensing gal axies in SDSS ([Shen et al.1 
I2003t IHvde fc Bernardill2009al) . This is due in part to the 
curvature in the r c -M* relation; SLACS is dominated by 
galaxies with M* > 10 1:l Mq whereas lO n M0 is the mid- 
point of the data fit in the SDSS samples. Indeed, we see 
in Figure [2] that the SLACS lenses follow the high-mass 
end of the quadratic fit of Hvdc & Bernardi ( 2009a) rea- 
sonably well. Additionally, SLACS is effectively a sample 
selected on velocity dispersion and the typical velocity 
dispersion at fixed effective radius is likely larger than 
non-lensing SDSS galaxies; this would lead to higher in- 
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TABLE 2 
M„ and a e / 2 — M» Relations 



Y 



X 



Slope a Intercept b Scatter 



e/2 
e/2 



M* 0.24 ±0.02 

M, 0.18 ±0.03 

M« 0.89 ± 0.04 

M» 0.81 ±0.05 



2.32 ±0.01 

2.34 ±0.01 0.04 ±0.01 

0.52 ±0.02 

0.53 ±0.02 0.05 ±0.02 



Note. — Fits arc of the form log Y 



log X + b 



with M* in units of 10 Mq, fr e / 2 in units of km s~ , 
and r e in units of kpc. Each fit is performed twice, ei- 
ther including intrinsic scatter or assuming zero intrinsic 
scatter. Note that the inclusion of intrinsic scatter has 
a significant affect on both the size-mass and velocity 
dispersion-mass relations. These fits are for a Chabricr 
IMF, but the slope is unaltered assuming a Salpctcr IMF. 

ferred stellar masses at fixed effective radi us (see Figure 
[p, as is seen in our data. Furthermore, iTortora et al.l 
(|2009t ) find a slope of 0.73 ± 0.12 for the r e -M* relation 
of massive (M» > 10 111 M Q local ETGs, completely con- 
sistent with our results. 

There are two components of the velocity dispersion 
selection function in the SLACS sample. The first comes 
from the lensing cross section which scales approximately 
with a 4 . The second one comes from the selection func- 
tion of the SDSS survey. Firstly, SDSS is a flux-limited 
sample so that high luminosity, and therefore high a, 
galaxies are over-repre sented because they are visible 
over a larger volume (Hyde fc Bernard! I2009af) . Sec- 
ondly, SDSS has finite resolution and thus ultra-compact 
galaxies are d ifficu lt to identify (e.g. JTruiillo et al.ll2009t 
iTavlor et al.ll2009t [Stockton et al.|[2010lh 

We can correct for the lensing bias to make the SLACS 
sample directly comparable to the parent SDSS sample 
by weighting each galaxy's contribution to the poste- 
rior distribution function by an exponent proportional 
to <7~ 4 (in \ 2 terms this would be equivalent to weight- 
ing each galaxy's contribution to the x by the same 
factor). Additionally, we can weight each galaxy by the 
volume in which it coul d be observed to p r ovide a more 
direct comparison with IHvde fc Bernardil (p0~0 9a) . We 
find that these weighting schemes alter the fits that do 
not include intrinsic scatter significantly more than the 
fits with intrinsic scatter; the fits with intrinsic scatter 
yield consistent results with and without weighting, and 
we therefore quote the unweighted fits throughout. 

3.2. Correlations between Stellar, Dynamical, and Total 

Mass 

Bivariate correlations between mass estimators are 
shown in Figure [3] and the parameters of linear fits to 
these relations are given in Table [3] The linear correla- 
tion between dimensional mass Md; m and lensing (or to- 
tal) mass M ro / 2 indicates that the virial coefficient is con- 
stant over the range in mass probed by the SLACS sam- 
ple, in agreement with our previous result (|Bolton et al.l 
2008b). The average value of the dimensionless parame- 
ter akin to the virial coefficient, log c C 2, defined by 



log M ro/2 = log 



c C 2r c cr i 



e/2 



2G 



is found to be 0.53 ± 0.09 and the scatter is 0.06 ± 0.01 
dex; both of these are consistent with our previous 
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Fig. 2. — The r c -M, relation for SLACS lenses. The solid black 
line is a linear fit to the relation; the dotted black lines indicate 
the intrinsic scatter and the gray band indicates the quadrature- 
sum of the scatter and the uncer tainty on the linea r fit. The blue 
dotted line is the relation fit by Shcn ct al. (2003) while the red 
dotted line is the linear fit from Hyde & Bcrnardi (2009a) and the 
red dashed line is their quadratic fit. 



measurement dBolton et al.|[20~08b[ ). The uniformity of 
the virial coefficient (i.e., the very small intrinsic scat- 
ter) does not imply exact scale-invariance of the mass- 
dynamica l structure o f early type galaxies. In fact, as 
shown by Nip oti et al.l ([20081 ). the observed virial coeffi- 
cient can be reproduced by a variety of two component 
mass models with a broad distribution of central dark 
matter fractions. The uniformity of the virial coefficient, 
however, restricts the possible range of acceptable mod- 
els for ETGs; in particular, models with extreme orbital 
anisotropics or that depart significantly fro m an isother- 
mal t otal mass density profile are ruled out ([Nipoti et al.l 
12001 . 

The correlation between total mass and stellar mass 
has the same amount of intrinsic scatter as the one with 
dimensional mass, once the larger errors associated with 
stellar mass are accounted for. However, the slope of the 
correlation differs significantly from unity. This is consis- 
tent with the well-known "tilted" slop e of the correlation 
between dy namical an d stellar mass (iFaber et al.l Il987t 
Cappcllari ct al. 2006; Gallazzi et al. 2006; Bundy et al. 
2007HRettura et al.ll2006HGraves & FaberlfeoiOD and can 



be explained in terms of a varying dark matter frac- 
tion and/or stellar IMF normalization with mass, in the 
sense that more massive galaxies have a higher fraction of 
(baryonic or otherwise) dark matter. Interestingly, how- 
ever, the small intrinsic scatter indicates that at fixed 
mass the dark matter fraction and/or stellar IMF nor- 
malization are tightly constrained. We will return to 
these points in § 13.41 

3.3. Correlations with the Slope of the Total Mass 
Density Profile 

We first examine the overall distribution of slopes of 
the mass density profile in a joint-framework wherein the 
7' of all of the early-type lenses are assumed to be drawn 
from a normal distribution parameterized by an average 



7 and a dispersion a ' . We find 
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Fig. 3. — Bivariate correlations between dimensional (Mdi m )i total (M r / 2 ) and stellar mass (M»). The best fitting linear relations 
are shown as dotted lines. Their coefficients are given in Table [3] Note the linear relation between dimensional and total mass, and the 
non-linearity of the other two relations. These are consistent with a constant virial coefficient and an increase with mass of the dark matter 
content and/or stellar initial mass function normalization. 

mining the power-law slope from a e / 2 and ctsie without 
n eeding to perform Jeans modeling. 

iKoopmans et al.l (J2009T ) found that the power-law slope 
did not correlate strongly with many of the global galaxy 
observables, including redshift, the ratio of the Einstein 
and effective radii, the central lensing mass, and ctsie- 
This updated analysis does not substantially change 
these results since we are only introducing 20% more 
systems. However, we now also investigate correla- 







TABLE 3 






Correlations Between Masses 




Y 


X 


Slope 


Intercept 


Scatter 


M dim 


M ro/2 


0.97 ±0.02 


0.51 ±0.03 




M dim 


Mr-e/2 


0.98 ± 0.03 


0.50 ±0.04 


0.06 ±0.01 


M* 


Mr-e/2 


0.81 ±0.03 


0.35 ±0.04 




M* 


Mr-e/2 


0.80 ± 0.04 


0.36 ±0.05 


0.03 ± 0.02 


M* 


Mdim 


0.80 ± 0.04 


-0.01 ±0.06 




M* 


Mdim 


0.79 ± 0.04 


0.01 ±0.08 


0.04 ±0.02 



Note. — Correlation between masses. All fits are done 
in logarithmic scales and using masses in units of 10 Mq, 
to reduce covariance. Fits without and with intrinsic scat- 
ter are given for each pairwisc combination. For example, the 
first line contains the results of fitting logMdi m /10 Mq — 
a logM 7 . c /2/IO Mq -t-b, where a and b are the slope and inter- 
cept, respectively. The second line adds an additional Gaus- 
sian component with average zero to represent intrinsic scatter. 
The width of the Gaussian intrinsic scatter is 0.06 ± 0.01 dcx 



and a' = 0.16 ± 0.02, consistent with and extending the 
results of IKoopmans et al.l (|2009D . The SLACS early- 
type lenses appear to be slightly super- isothermal (i.e., 
the density profiles are typically steeper than isother- 
mal by ~ 5%) with an intrinsic spread of ~ 10%. We 
note, however, that we have assum ed that th e ga laxies 
are isotropic; as discussed in Koopmans et al. (2009), a 
modest amou nt of radial an isotropy (f3 < 0.5, consis- 
tent with iGer hard et aT. 2001) is sufficient to produce an 
isothermal slope, 7 = 2. We ca nnot directly probe the 
anisotropy with our data (but see IKoopmans et al.ll2009L 
which provides an indirect constraint) and we therefore 
impose the isotropic model. 

It was previously found that the ratio of the ob- 
served stellar velocity dispersion and the SIE model 
velocity dispersion, /ste = cr e /2/ (T sm, is strongly cor- 
related with the logarithmic density slope (|Treu et al.l 
2009). Our updated analysis confirms this result (Fig- 
ure |1|) and the best-fit linear relation, taking into ac- 
count the covariance between /sie and 7', is given by 
i - 2 = (2.67 ± 0.15)(/ S ie - 1) + (0.20 ± 0.01). Note 
that this trend is a natural consequence of how 7' is de- 
termined and does not provide any significant physical 
insights; instead, it provides a useful shortcut for deter- 



tions with r c , <r e / 2 , and the central surface mass den- 
sity Etot = M^^Ae- We find non-negligible correlations 
(non-zero slopes with greater than 3-ct significance) with 
r c and E to t but no clear trend with a e / 2 (Figure [5]); the 
correlation with E to t is the tightest and most significant 
(Table E]). This is expected, since a steeper mass den- 
sity profile implies a higher central surface mass density, 
and explains at least in part the intrinsic dispersion in 
the average mass density profile. However, the residual 
intrinsic dispersion (0.12, i.e. 6% in slope) indicates that 
there may be additional observables that correlate with 
the inferred slope. One such parameter could be local 
environment, due to tidal effects on the outer halos or 
to contamination to the lensing convergence by external 
mass along the line of sight. The former has been ten- 
tatively detected usi ng the SLACS sample at ma rginal 
levels of significance ([Auger 200cJ iTreu et aT1l2009f), and 
it ha s been observed in clusters (e.g.. iNatarajan et al.1 
2009). The latter does not seem to be significant in the 
SLACS sample, as inferred from the minimal level of mis- 
alignment between the major axis of the light and mass. 
Another element that may contribute to the scatter is 
anisotropy of the stellar orbits. However, higher preci- 
sion measurements will be required to determine whether 
the residual scatter is stochastic in nature and/or if there 
are residual and undetected small systematic trends. 

There is tentative evidence for a slight anti-correlation 
between 7' and the total and stellar masses, as one might 
infer from the trends with radius and surface mass den- 
sity, although neither of these anti-correlations between 
slope and mass are statistically significant given our sam- 
ple size and data quality. Nevertheless we note that an 
anti-correlation between m ass and density slope can aris e 
as a result of mergers (e.g.. lBovlan-Kolchin fc Ma 2004), 
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Fig. 4. — The relationship between the ratio of the observed 
stellar velocity dispersion and the SIE model velocity dispersion, 
/sie = °"e/2/°"SlEi an d the logarithmic density slope 7'. The 
correlation is very tight, with evidence for little intrinsic scatter 
(o'int = 0.02). Note that although the errors on / and 7' are 
strongly correlated our analysis takes this correlation into account 
and finds a significant excess correlation. The tightness of the 
relation implies that /sie can be a useful proxy for 7' without 
performing the joint lensing and dynamics modeling. 
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Fig. 5. — Correlations between r c , cr e /2, and Etot (see Table [4] 
for the linear fits). The correlation with o- e / 2 is not significant al- 
though the r e and Etot correlations are found at greater than 3-cr 
significance. All correlations show significant scatter, and the ef- 
fective radius correlation is somewhat weaker than the central sur- 
face mass density correlation. This latter correlation is expected; 
steeper power-law slopes imply higher central concentrations, and 
higher central concentrations imply increased central surface mass 
densities. 



although analytic models of halo collapse may predict a 
positive correlation between the m ass and central density 
slope of the dark matter halo (e.g.. lDel Popolo fc Kroupal 
I20M) . 

We have previously combined the SLACS lenses with 
a higher-redshift set of l enses from the Lensing Struc- 
ture and Dynamics fLSD: rTreu fc Koopmansll2004D sam- 
ple and found marginal evidence that the central mass- 
density slope evolves with redshift (jKoopmans et all 



TABLE 4 
Correlations with 7' 



X 



Slope 



Intercept 



Scatter 



-0.41 d= 0.12 0.39 ±0.10 0.14 ±0.02 
0.07 ±0.08 -0.12 ±0.21 0.17 ±0.02 
0.85 ±0.19 -0.47 ±0.12 0.12 ±0.02 



CT e/2 
Stot 



Note. — r c is in units of kpc, o- e / 2 is in units of 
100 km s" 1 , and E t ot is in units of 10 11 M kpc -2 . 

2006). However, the relationship between 7' and E to t 
was not explicitly accounted for in that analysis; a full 
investigation of the evolution of 7' including this effect 
will require an expanded sample of high-redshift lenses. 

3.4. Central mass-to-light ratios and dark matter 
fraction correlations 

The SLACS dataset presents a unique opportunity to 
investigate the central mass-to-light ratio and dark mat- 
ter fraction in galaxies beyond the local universe. We 
first look at the relationship between M/L and six other 
parameters: B-band luminosity at z = 0, V-band lu- 
minosity at z = 0, cr e / 2 , Mdim, M*, and M r e /2- These 
trends are shown in Figure [5] and listed in Table [5] It is 
clear that there is a significant trend between all of these 
parameters and total M/L, while M*/L is, at the 2-cr 
level, independent of the quantities investigated, and we 
also note that M*/L and the total M/L correlate more 
strongly with <T e / 2 than with mass or l uminosity. These 
result s are i n excellent agr e emen t with iCappellari et all 
(|2006f) and iTortora et all (|2009f) . who investigate the 
stellar M/L and dynamical M/L for samples of E and 
SO galaxies, a nd with results from the FP and M»P of 
SPSS galaxies (Tfvde fc Bernardii 2009b; Graves fc Faberl 
120101 ). However JGrillo fc Gobatr^OlOl ) find a somewhat 
steeper trend between M* and M*/Lb for the SLACS 
lenses (a — 0.18), which we attribute to differences in the 
stellar mass determin ation due to assum ptions about age 
and metallicity (e.g.. lAuger et al1l2009h . For very mas- 
sive ETGs like the ones in the SLACS sample, the differ- 
ences in stellar population properties are not sufficiently 
large to account for large changes in the stellar M/L with 
a fixed IMF. Conversely, the total mass-to-light ratio 
clearly increases with total stellar mass, possibly as a re- 
sult of increased dark matter or varying stellar IMF. The 
stellar M/L relations are all consistent with no intrinsic 
scatter, indicating a remarkable homogeneity in the stel- 
lar populations of these galaxies. The total M/L trends, 
on the other hand, exhibit 0.07-0.09 dex of intrinsic scat- 
ter; this is expected since the various parameter planes 
(§ HJ demonstrate that three pa rameters are required t o 
adequately describe ETGs (e.g. JGraves fc Faberl 1201 Of ). 

In Figure [7] we show the trends in the projected dark 
matter fraction with the total mass inferred from lensing, 
the stellar mass, the observed stellar velocity dispersion 
from SDSS, the effective radius, and the B- and V-band 
luminosities evolved to z — 0. Linear fits to these re- 
lations are provided in Table El The most significant 
trends are with effectiv e radius and total mass (also see 
iNapolitano et al]|201Q[ ). indicating that these parameters 
govern the central dark matter fraction. These parame- 
ters are also consistent with having no intrinsic scatter 
at the 95% confidence level (i.e., a- m t is within 2-er of 
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TABLE 5 
M/L ] Linear Relations for SLACS Lenses 







Stellar M/L 






Total M/L 




X 


a 


b 


Omt 


a 


b 


Omt 


L B (B-band) 


0.02 ± 0.04 


0.59 ±0.03 


0.02 ±0.01 


0.19 ±0.06 


0.85 ± 0.05 


0.09 ±0.01 


L B (V-band) 


0.01 ± 0.04 


0.49 ±0.03 


0.01 ±0.01 


0.19 ±0.06 


0.75 ±0.05 


0.09 ±0.01 


L v (B-band) 


0.03 ± 0.04 


0.58 ±0.04 


0.02 ±0.01 


0.21 ±0.06 


0.82 ± 0.05 


0.09 ±0.01 


L v (V-band) 


0.02 ± 0.03 


0.49 ± 0.03 


0.01 ±0.01 


0.20 ± 0.06 


0.71 ±0.05 


0.09 ±0.01 


CT e/2 (B-band) 


0.26 ±0.16 


0.50 ±0.07 


0.02 ± 0.01 


0.89 ± 0.23 


0.63 ± 0.09 


0.09 ±0.01 


°e/2 (V-band) 


0.19 ±0.14 


0.42 ± 0.06 


0.01 ±0.01 


0.81 ±0.23 


0.56 ±0.09 


0.09 ±0.01 


M dim (B-band) 


0.03 ± 0.03 


0.59 ±0.02 


0.02 ± 0.01 


0.24 ± 0.05 


0.83 ± 0.03 


0.08 ±0.01 


M dim (V-band) 


0.02 ± 0.03 


0.49 ± 0.02 


0.01 ±0.01 


0.23 ± 0.04 


0.72 ± 0.03 


0.08 ±0.01 


M, (B-band) 


0.07 ±0.04 


0.58 ±0.02 


0.02 ±0.01 


0.22 ± 0.06 


0.91 ±0.02 


0.09 ±0.01 


M, (V-band) 


0.05 ± 0.03 


0.48 ±0.01 


0.01 ±0.01 


0.21 ±0.06 


0.81 ± 0.02 


0.09 ±0.01 


M r e /2 (B-band) 


0.04 ± 0.03 


0.60 ±0.01 


0.02 ±0.01 


0.27 ±0.04 


0.93 ±0.01 


0.07 ±0.01 


M r„/2 (V-band) 


0.03 ± 0.03 


0.50 ±0.01 


0.01 ±0.01 


0.26 ± 0.04 


0.83 ±0.01 


0.07 ±0.01 



NOTE. — Fits are ot the form log[ M/L / (M/L)q] = a*log[Jf] + b and the stellar M/L is determined 
assuming a Chabricr IMF. Lb and Lv are in units of 10 Lq, <T e / 2 ls in units of 100 km s — , and the masses 
(M„, M dim , and M rc/2 ) are in units of 10 11 M . 



0), although the errors are large. We note that part 
of this scatter may also be due to the inadequacy of a 
linear fit; this is most clear for the Salpeter relations, 
which would indicate a physically impossible negative 
dark matter fraction if the linear trend is extrapolated 
below a stellar mass of ~ 10 105 M Q . A more appropriate 
model will be examined in a future paper (Auger et al. 
2010, in preparation). 

4. CORRELATIONS IN THREE OR MORE DIMENSIONS 

The SLACS lens galaxies have previ ously been 
show n to lie on t h e Fun damental Plane (jTreu et al.l 
120061: IBolton et all [2008b') and on the Mass Plane 
(jBolton et al.l 120071 12008b). We now revisit these 
relations with the enlarged SLACS sample and self- 
consistently evolved luminosities, superceding our previ- 
ous analysis. We also investigate other scaling relations 
involving stellar mass, such as the stellar mass including 
the Stellar Mass Plane (M*P) and a new correlation in 
a higher dimensional parameter space, which we call the 
Fundamental Hyper-Plane. 

4.1. Fundamental and Mass Planes 
We fit the parameter plane relations with the form 

log r e = a pp log a e/2 + /? pp log A + 7 PP , (2) 

where A represents the average surface brightness within 
r e , the average stellar mass surface density within r c , or 
the average total mass surface density within r c /2. The 
units used for the fits are: kpc for r c , 100 km s^for a e / 2 , 
10 9 Lq for the luminosity, 10 9 M for the stellar mass, 
and 10 10 Mq for the total mass. The intrinsic scatter is 
given in units of log r c . The inferred parameter planes 
are shown in Figure [5] and illustrate the small intrinsic 
scatter found in these relations. 

We find that the FP relation is s omewhat tilted with 
respect to previous analyses fe.g.. IBolton et all 12008b) 
if we allow the intrinsic scatter to be a free parameter 
of our fit (see Table [7]). If we do not fit for intrinsic 
scatter (that is, if we impose that the intrinsic scatter 
is zero) we recover a FP con s istent with other determ i- 
nations (jBolton et al.l l2008bl: iHvde fc Bernardll2009bO . 
This illustrates the importance of explicitly accounting 



for intrinsic scatter in the fit, even when it is small like 
in the case of the FP. As a further consistency check 
we repeated the fit of the FP parameters applying the 
<r~ 4 scaling to account for the velocity dispersion selec- 
tion of the SLACS sample, and we find that the changes 
are insignificant. This is consistent with the findings of 
IHvde fc Bernardil (|2009bl ). since our fitting method is 
closer to their "direct" fit than to their orthogonal fit. 
The non-zero scatter of the FP confirms previous results 
and is consistent with being due to a combination of stel- 
lar population differences and structural differences. The 
availability of stellar mass and total mass diagnostics al- 
lows us to break this degeneracy as we discuss in the rest 
of this section. 

The coefficients for the M*P plane are independent of 
the choice of IMF for the SPS models (the normaliza- 
tion term changes by « 0.25, as is expected for Chabrier 
and Salpeter IMFs when /3 m — 1). It has previously 
been found that the M*P (where M* is inferred from SPS 
models) lies closer to the virial plane (a = 2; /3 = — 1) 
than the FP (|Hvde fc Berna rdi 2009b]), although we find 
in our data that the FP and M^P are approximately 
aligned. We note, however, that the M*P is consistent 
with having no intrinsic scatter (<7i, lt = 0.020 ± 0.014) 
while the FP has intrinsic scatter of <T; nt — 0.049 ± 0.009 
dcx; this implies that the scatter in the FP is largely 
driven by small differences in stellar populations (e.g., 
age or metallicity) for massive galaxies, rather than from 
structural properties (i.e., differences in dark matter con- 
tent or the virial coefficient at a fixed size and velocity 
dispersion). 

The MP is also consistent with having no scatter (Ta- 
ble [7]) but is substantially mis-aligned with the FP and 
M*P; instead, the MP is approximately aligned with the 
virial p lane (as was found in Bolto n et al.l 1200 7, 2008b; 
also see lKoopmans et al.ll2009l where this plane included 
7'). The offset with respect to the FP and M*P is consis- 
tent with the bivariate relations (e.g., the middle panel 
of Figure [3] and the total M/L trends of Figure [6]), while 
the slight offset from the virial plane may be related to 
the systematic deviation of the central mass density pro- 
file from isothermality (e.g., Figure[S]) or anisotropy. The 
MP is also found to be consistent with no intrinsic scat- 
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TABLE 6 
/dm Linear Relations for SLACS Lenses 







Chabrier IMF 






Salpeter IMF 




X 


a 


b 


Tint 


a 


b 


0"int 


Lb 


0.16 ±0.06 


0.47 ±0.05 


0.05 ±0.02 


0.28 ±0.10 


0.08 ±0.08 


0.09 ±0.03 


Lr 


0.16 ±0.06 


0.46 ±0.05 


0.05 ±0.02 


0.28 ±0.10 


0.04 ± 0.09 


0.09 ± 0.03 


°"e/2 


0.46 ±0.22 


0.40 ± 0.09 


0.06 ±0.02 


0.80 ± 0.44 


-0.05 ±0.18 


0.11 ±0.03 


r e 


0.28 ±0.05 


0.36 ±0.05 


0.03 ±0.02 


0.49 ±0.10 


-0.13 ±0.09 


0.06 ± 0.03 


M, 


0.13 ±0.06 


0.54 ±0.03 


0.06 ±0.02 


0.23 ±0.11 


0.14 ±0.07 


0.10 ±0.03 


M r c /2 


0.20 ±0.04 


0.54 ±0.02 


0.03 ±0.02 


0.36 ± 0.07 


0.20 ±0.03 


0.06 ± 0.03 



Note. — Fits are of the form /dm — a*log[X] + b. Lb and Ly arc in units of 10 10 Lq, 
in units of 100 km s~ , r c is in kpc, and M* and M r / 2 arG m units of 10 Mq. 



ter, although the errors are large. 



4.2. The Fundamental Hyper-Plane 

The M*P and MP are found to have little intrinsic 
scatter but very different orientations in size-velocity 
dispersion-mass space. We would like to know whether 
the remaining scatter in these planes is due to a non 
trivial relationship between M ro / 2 and M*. We therefore 
explicitly explore the relative importance of the stellar 
and total mass on the structure of ETGs. In particular, 
we consider the relationship between the effective radius, 
velocity dispersion, central stellar mass, and central total 
mass by assuming a relation of the form 

log r c = a hp log a e/2 + /3^log M rc/2 + 7 h Plog M* + S h ? 

(3) 
where M* is now the stellar mass within half the effective 
radius to ensure consistency with the lensing determined 
total mass. An initial fit to this relation finds that a small 
number of galaxies, six in total, are significant outliers 
(greater than 3-(7i n t, where o~i n t = 0.03 for the initial 
fit). We re- fit the relation without these objects and 
find a substantially tighter relation with approximately 
the same coefficients but substantially decreased scatter. 
We therefore suspect that these objects have aberrant 
velocity dispersions and have therefore excluded them 
from our analysis, as discussed in Section 2. 

The best-fit relation for the Fundamental Hyper-Plane 
(FPH), with r c measured in kpc, a e /2 in 100 km s , and 
both stellar and total mass in 1O 1O M , is given by 

log r c = (-0.91 ± 0.10)log <r e/2 + (0.69 ± 0.04)log M rc/2 + 

(0.11 ± 0.06)log M* + (0.23 ± 0.03) (4) 

with intrinsic scatter a ln t = 0.007 ± 0.005, and the fit is 
shown in Figure [5] The dominant terms are the central 
velocity dispersion and central total mass, while the stel- 
lar mass only has a marginal role and is consistent with 
being unimportant at the 2-er level (that is, the coefficient 
for the M* term is within two-sigma of zero) . 

5. DISCUSSION 

Before discussing our results in the broader context of 
studies of the structure formation and evolution of ETGs, 
it is essential to test whether the SLACS ETGs are in- 
deed representative of the overall population. In previous 
papers we have shown that SLACS lens galaxies cannot 



be distinguished from control samples of ETGs selected 
from the SDSS archive t o have the same s t ellar velocity 
dispersion and redshift (JTreu et ail 2006; B olton et al.l 
l2008al lAuger et al.ll2009T ) . In this paper we further refine 
these tests by considering the complete SLACS sample of 
73 ETGs, and testing a number of correlations, includ- 
ing the FP, the M*P and their projections, the M*-r e and 
M*-ct correlations. Once again we do not find any sub- 
stantial difference between the correlations inferred for 
our sample and for SDSS-selected samples of non-lenses. 
Having found no significant evidence of a systematic dif- 
ference, we can safely assert that SLACS lenses are rep- 
resentative of the entire population of massive ETGs and 
proceed to interpret our results. 

5.1. The "bulge-halo" conspiracy, and implications of 
(small) departures from it 

The first key result of this study, building on previ- 
ous SLACS papers, is the precise measurement of the 
so-called "bulge-halo conspiracy" and its tightness. In 
short, although mass clearly does not follow light and an 
extended dark matter halo is needed to reproduce simul- 
taneously the lensing and dynamical constraints, the two 
components add up to form almost exactly an isothermal 
total mass density profile in the inner regions of galaxies 



(fh 



with 7' = 2). This is remarkable, since 



neither component is a single power law, and there is 
no simple fundamental reason why this should be the 
case, although general dynamical arguments based on 
incomplete violent relaxation suggest that the is other 
mal sphe re is a form of dynamical attracto r (e.g.. | Gun: 
1977; Dckcl et al. 1981; van Albada 1982; Bertschingc: 
ll985J ILoeb k Peebles! 120031 ). In addition, models based 
on simple prescript ions for baryonic condensation (e.g., 
IGnedin eT al. 2004) seem to suggest that it is possible 
to obtain close-to-isothermal total mass density profiles 
star ting from cosmo logically motivated dark matter ha - 
los (| Jiang fc Kochanek|[2007l: Humphrey fc Buoteil2009ft . 
The resulting isothermal profile cannot be explained 
based purely on dissipationless processes. In fact, dis- 
sipationless processes in a cosmological setting would 
tend to produc e inner density pro files close to 7' = 1 
or even flatter (jNavarro et al.ll2~004D . However, once the 
isothermal p rofile is established via dissipational pro - 
cesses (e.g., iCiotti et al.l 120071 : iRobertson et al.l [2006). 
collisionless "dry" mergers preserve it quite accurately, 
introducing a small amount of scatter , consistent with 
the observed value (|Nipoti et al.ll200"9alh An interpreta- 
tion of the bulge-halo conspiracy and its possible origins 
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Plane 



TABLE 7 
Parameter Planes for SLACS Lenses 



, pi> 



P 



)'!' 



7" 



FP 1.189 ±0.141 -0.885 ±0.041 

M,P 1.191 ±0.221 -0.971 ±0.073 
MP 1.829 ±0.133 -1.301 ±0.061 



FP 1.020 ±0.203 

M„P 1.185 ±0.214 
MP 1.857 ±0.136 



-0.872 ± 0.052 
-0.952 ± 0.074 
-1.279 ±0.065 



-0.185 ± 0.047 

0.257 ±0.088 

-0.301 ± 0.055 

-0.108 ± 0.076 

0.261 ± 0.085 

-0.312 ±0.056 



0.049 ± 0.009 
0.020 ±0.014 
0.013 ±0.010 



Note. — Fits are of the form given in Equation[2] with r c in units of kpc, a e i 2 
in units of 100 km s~ , V-band luminosity and stellar mass in units of 10 Lq 
and 10 Mq respectively, and total mass in units of 10 Mq. The first three fits 
are without i ntrinsic scatt er while the latter three explicitly include scatter using 
the model of lKellvl 1 120071) . 



are discussed at length s in previous SLACS papers (e.g., 
iKoopmans et al.ll2006[ ) and will not be repeated here. 

However, the prec ision achieved in this and p revious 
SLACS studies (e.g.. IKoopmans et al.ll2006l 120091 ) allows 
us to highlight the importance of the small but signifi- 
cant observed departures from the bulge-halo conspir- 
acy. Firstly, the average mass density profile is not ex- 
actly isothermal but slightly steeper, 7' = 2.078 ± 0.027 
(this is marginally steep er than the slopes found by 
iHumphrev fe B uote 20091 who use X-ray temperature 
and density profiles to find the best-fit power-slope for 
the mass distribution of four ETGs). Secondly, there is 
evidence for non-negligible intrinsic scatter 0.16 ± 0.02. 
Third, there is tentative evidence for a mild dependency 
between 7' on galaxy properties, such as radius or central 
mass density. 

The first fact has important implications for gravi- 
tational lens studies, particularly those trying to infer 
cosmography from gravitational time delays. Given the 
known degeneracy between slope of the m ass density 
profile and time delays (e.g., Wucknitz 2002), if one as- 
sumes an isothermal prior then the inferred Hubble Con- 
stant will typically be biased low by 10% if no ot her di- 
rect measure of the mass slope is available (see lOgurl 
[20071 : IDobke&Kinl 120061 : IKochanekl [20061 : iTreul I2010L 
and references therein). Additionally, due to the in- 
trinsic scatter, estimates of the Hubble Constant from 
a single lens can be off by 20% if an isothermal model 
is assumed and not independently constrained. Note, 
however, that these arguments are based on our assump- 
tion of isotropic orbits; anisotropy could chang e these bi- 
ases considerably (e.g.. IKoopmans et al.ll2009T i and must 
be constrained in more detail in future work. Neverthe- 
less, additional external information on the mass density 
slope, such as t hat inferred from multiply-imaged ex- 
tended sources ([Warren fe Dvd 120031: ISuvu et al. 120061) 
and from stellar kine matics ( Treu fc Koopmana 120021 : 
IKoopmans et al.l l2003f ). can help to break these degen- 
eracies and pr ovide robust estim ates of cosmological pa- 
rameters (e.g. JSuvu et al.|[20 09). 

From the point of view of galaxy formation the intrin- 
sic scatter and the correlations between 7' and galaxy 
pro perties are the most inter esting new elements (also 
see [Humphrey fc Buotd 120091 ) . The tightness imposes 
a constraint on the number of major merging events, 
as well as on the diversity of formation histories. We 
will return to this point later in § 15.31 in the context 



of the tightness of the empirical correlation. The de- 
pendency of mass density slope on galaxy parameters 
is another piece of evidence supporting the idea that 
ETGs are not a homologous family from a structural 
point of view. Both facts provide an interesting con- 
straint for numerical simulations that have sufficient res- 
olution and bar yonic physics to s imulate the inner r egions 
of ETGs (e.g.. [Robertson et al.l 120061: iNaab et aLl 12007 ' 



aDJlQOl 

BrT T2(r7nF . 



IGonzalez-Garcfa et al.l l2009t lLackner fc Qstrikerll201 

Small departures from regularity in the end may be 
the key to understandi ng the details of ETGs formation 
(jKormendv et al.ll2009T ). 

5.2. Non-triviality of empirical correlations 

The second key property of ETGs that we quantify in 
this study is the degree to which their internal structure 
changes as a function of mass or size. There are many 
ways to express this, including the observed "tilt" of the 
FP and the mass-dependent correlation between stellar 
mass and dynamical mass. The lensing observables from 
the SLACS sample allow us to add additional informa- 
tion and breaking some of the degeneracies in the inter- 
pretation of these trends. We find that, over the range of 
masses probed, the total amount of mass is a non-linear 
function of stellar mass, while the virial coefficient is ap- 
proximately constant. Similarly, we find that the trends 
cannot be due to changes in stellar mass-to-light ratio for 
a fixed IMF. In fact, the M*/L inferred from SPS models 
is approximately th e same for all SLAC S ETGs (see also 
iTortora et all 120091 : iGrillo et al.l 120091 ) . consistent with 
the homogeneity of old stellar populat ions present in 
massive ETGs (e.g.. lThomas et al.ll2005l) . The fact that 
the total M/L varies strongly as a function of mass, veloc- 
ity dispersion, and luminosity means that more massive 
galaxies have higher total M/L than less massive sys- 
tems. This suggests that either the central dark matter 
fraction or the IMF is a stro ng fun ction of mass (e.g., 
iTreu et al.H2010HTortora et al.ll2009h . In other words, it 
is the central fraction of dark matter, either baryonic or 
non-baryonic, that increases with galactic stellar mass. 

These trends have long been connected with the pres- 
ence of some characteristic scale in the formation pro- 
cess of ETGs, that may be due to baryonic physics. For 
example, the increase in cooling timescales for the hot 
gas going from the lower mass to the higher mass ETGs 
significantly changes the efficiency of converting baryons 
into stars and therefore affects the dark matter fraction in 
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Fig. 6. — Relations between the stellar (crosses) and total (points) M/L in the B- (blue) and V- (red) bands for the SLACS lenses. A 
Chabrier IMF has been assumed for the stellar M/L. The dotted lines are linear fits to the stellar M/L relations and dash-dotted lines are 
linear fits to the total M/L relations. 
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the central regions fe.g. JRobertson et al.ll2006l and refer- 
ences therein). Additionally, changes in the modes of star 
formation that are responsible for variations in chemical 
abundances coul d perhaps induce vari ations in the stel- 
lar IMF (but see iGraves fc Faberll2010t for a discussion) . 
However, collisionless processes can also contribute to 
the emergence of changes in structural properties with 
mass. For example, the fraction of dark matter within 
the cylinder of radius equal to a fixed fraction of the ef- 
fective radius changes during dry merge rs as a resu lt of 
the increase in effective radius (Ni poti et alJl2009aj ). 

Interestingly, the observed trends cannot be explained 
by a single phenomenon (for example, an increase in star 
formation efficiency with total mass) as we know that 
ETGs occupy at least a two dimensional subset of pa- 
rameter space even when stellar population effects are 
excluded by purely structural correlations like the Mass 
Plane and (assuming robust SPS models and ignoring the 
unknown normalization due to the IMF) the Stellar Mass 
Plane. It appears that velocity dispersion as well as size 
(or dynamical mass) are needed to fully specify the dy- 
namical properties of an ETG. Furthermore, neither the 
MP nor the M*P appear to have intrinsic scatter. Thus, 
remarkably, it appears that two parameters are not only 
necessary but also sufficient to fully specify the internal 
structure of a massive ETG within our observational er- 
rors. Finally, when we construct a direct relationship 
between the structural parameters (effective radius, ve- 
locity dispersion, stellar mass, and central total mass) of 
the SLACS lenses (EquationUJ) we find that the relation- 
ship is nearly independent of the stellar mass. In terms 
of observables for non-lens galaxies, the driving param- 
eters are size and stellar velocity dispersion, not stellar 
mass. 

5.3. Tightness of empirical correlations 

The third key feature of ETGs addressed in this study 
is the tightness of the empirical correlations. Traditional 
non-lensing correlations, such as the FP, have sm all but 
non-zero intrinsic scatt er (e.g.. | Jorgensen et al.l 119961 : 
IHvde fc Bernardll20n9bl ). IGraves fc Faberl (12010T ) have 
recently emphasized the presence of intrinsic scatter in 
the M*P relation at the level of 0.02-0.03 dex, discussing 
several interpretations in terms of diversity in stellar IMF 
and dark matter content. This low level of intrinsic scat- 
ter cannot be ruled out by our data, although it should 
be noted that our sample is restricted to the most mas- 
sive systems and scatter may be mass-dependent. Our 
study of correlations involving total mass, including the 
MP, adds an important piece of evidence because the 
MP is independent of SPS stellar mass estimates and 
therefore can be used to break the degeneracy between 
the IMF and dark matter. Our study shows that the 
MP is consistent with having no intrinsic scatter, within 
the measurement errors of a few hundredths of a dex 
(omt = 0.013 ± 0.010). Future studies comparing the 
M*P with the MP to even higher degrees of precision 
will be able to quantify the contribution of IMF varia- 
tions to the intrinsic scatter of the M*P. 

The tightness of scaling relations is especially remark- 
able in a sce nario where evolution i s driven by major 
mergers (e.g. . Ivan der Wei et al.ll2009h . For example, dry 
mergers tend to move galaxies within the FP and MP 
correlations and preserve their tightness. However, dry 



mergers do not, in general, retain the tightn ess of the bi- 
varia t e projections of the paramet e r planes (INipoti et al.l 
12001 IBovlan-Kolchin et al.l [20061: INipoti et al.l I2009bfl . 
The properties of the progenitors must be finely tuned 
with the orbital parameter of the merger in order to pro- 
duce th e tight observed scaling relations. INipoti et al.1 
( 2009b) used these relations to show that only half of the 
mass in ETGs at z = can result from dry merging and 
dry merging ca nnot cause sup e r-mas s ive galaxies at high 
redsh ifts (e.g., iTrujillo et al.l 120061 : Ivan Dokkum et afl 
l2008h to evolve into present-day ETGs. 

6. SUMMARY 

We briefly summarize the most significant conclu- 
sions from our analysis of the early-type lenses from the 
SLACS survey. 

• The SLACS sample obeys all the standard correla- 
tions found for non-lensing ETGs, consistent with 
the hypothesis that it is representative of velocity 
dispersion selected ETGs. 

• Stellar kinematics and lensing data constrain the 
the slope of the total mass density profile (ptot oc 
r~ 7 ). The average slope is found to be close to, 
but slightly steeper than, isothermal, with (7') = 
2.078 ±0.027 and an intrinsic scatter of 0.16 ±0.02. 

• The total mass density slope 7' correlates with ef- 
fective radius (r e ) and central mass density (E to t) 
in the sense that denser galaxies have steeper pro- 
files. The residual intrinsic scatter is reduced but 
still significant (0.12 ± 0.02 for £ tot ). 

• Tight correlations are found between dimensional 
mass Mdim = 5<T G r ° , stellar (M„) and total mass 
(M tot ). The relationship between total mass and 
dimensional mass is found to be consistent with 
linear with very little scatter, implying that the 
virial coefficient of ETGs is constant over this mass 
range. The correlation between total (dynamical) 
mass and stellar mass is non linear (M* oc M°y 2 ), 
consistent with the hypothesis that the central 
CDM content and/or the normalization of the stel- 
lar IMF changes with mass. 

• Assuming a universal initial mass function (IMF) 
the stellar mass-to-light ratio is nearly constant 
over the range in masses probed by the SLACS 
ETGs. In contrast, the total mass-to-light ratio 
correlates strongly with lens properties; the most 
significant correlations are with the central veloc- 
ity dispersion and central total mass. As a result, 
the dark matter fraction within r e /2 is a monotoni- 
cally increasing function of galaxy mass and size. If 
the universal IMF assumption is relaxed, the trend 
could be explained at least in part by an increasing 
IMF normalization with galaxy mass. 

• The Mass Plane (MP), obtained by replacing sur- 
face brightness with surface mass density in the FP, 
is found to be tighter and closer to the virial rela- 
tion than the FP and the M*P, indicating that the 
scatter of those relations is dominated by stellar 
populations effects. 
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Fig. 7. — Relations between the projected dark matter fraction within half of the effective radius and hg, Ly, c e /2, r e , M„, and M r<! /2- 
Red points are for a Salpeter IMF and black points are for a Chabrier IMF; the most significant trends are with r B and M r / 2 . 
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Fig. 8. — FP, M„P, and MP relations for the SLACS lenses. All three planes show little intrinsic scatter ((Tint $ 0.05) and the FP and 
M*P are well-aligned. The MP is found to be offset from the FP and M*P but is approximately aligned with the virial plane. The red 
points indicate the 6 galaxies which are outliers on the Fundamental Hyper-Plane and are not included in the fits. 

a function of halo mass and as the result of dry and wet 
mergers, reproducing quantitatively the observed corre- 
lations and their tightness may be a significant challenge 
for galaxy formation models. A detailed modeling ef- 
fort will be presented in a follow-up paper, where weak- 
lensing data will be combined with the present data to 
strengthen the connection between the central part of the 
galaxies and the virial mass of the halos in which their a 
embedded (Auger et al. 2010, in preparation). 
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• We construct the Fundamental Hyper-Plane by 
adding stellar masses to the MP and find that the 
stellar mass coefficient is consistent with zero and 
there is effectively no residual intrinsic scatter. 

Our results demonstrate that the dynamical structure 
of massive ETGs is not scale invariant and that it is fully 
specified by size, stellar velocity dispersion, and total 
mass. Although the basic trends can be explained qual- 
itatively in terms of varying star formation efficiency as 
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